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Abstract 



We construct complete sets of invariant quantities that are integrals 
^ ' of motion for two Hamiltonian systems obtained through a reduction pro- 

I cedure, thus proving that these systems are maximally superintegrable. 

0^ ■ We also discuss the reduction method used in this article and its possible 

generalization to other maximally superintegrable systems. 

O ■ 1 Introduction 

' Superintegrable Hamiltonian systems, both in classical and in quantum mechan- 

ics, have attracted a lot of attention in the last decades, due to their physically 
relevant properties. The existence of closed periodic bounded orbits for the clas- 
^ . , sical models, and of accidental degeneracy of the energy levels for the quantum 

■ ones (removed by the existence of complete sets of commuting observables) , 

are among the most remarkable consequences of the rich symmetry structure 
possessed by these systems. 

Recent investigations lead to the construction of new superintegrable poten- 
tials, in flat spaces and curved spaces. As is well known, a maximally super- 
integrable system in a TV-dimensional configuration space has 2N — 1 integrals 
of motion, although only sets of N integrals (chosen in several different ways 
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among the 2N — 1 ones) are in involution. Among the classical examples of this 
type of Hamiltonians are well known systems: the Coulomb potential and the 
harmonic oscillator, whose dynamics is described by Bertrand's theorem P] (for 
a review see, for instance, [T7] ) 

Superintegrable systems are rare and not easy to construct. Nevertheless, the 
approach we describe here offers a systematic way to produce superintegrable 
systems. 

In this contribution, we are primarily interested in the reduction of super- 
integrable systems by using an adaptation of the Marsden-Weinstein reduction 
scheme |12j and in the study of the properties of the corresponding reduced 
systems. We will not analyze here in full generality and in all details the setting 
of this method. Instead, we will focus on two particular, but very interesting, 
cases. Reduction methods have been used as a common tool to derive new dy- 
namic Hamiltonians starting from free systems, defined in the original manifold 
(see, for instance, [3]). Instead, in our approach, we start with a dynamics 
that is not free, and study its reductions. The original potentials (Coulomb 
and harmonic oscillator) are already superintegrable, and their invariants are 
known. 

The first system, i.e. the Coulomb potential plus centrifugal terms, was 
analyzed in ,l9j by using an interesting ad hoc approach. We will present here 
a different point of view, which can be easily generalized to other models. We 
apply our approach to this model and to that one obtained as a reduction of 
the classical anisotropic harmonic oscillator. We recall that these systems, are 
also exactly solvable, i.e. their quantum spectrum can be computed by purely 
algebraic means [IH] . This paper is based on the work [15, . For a treatment of 
the quantum case, see 0. 

2 A geometric approach to the construction of 
superintegrable systems 

In this section, we briefly review the geometric settings motivating our results. 
Wc work on a 2iV— dimensional symplectic manifold (M, w). Recall that a mo- 
mentum map for the action of a group G on M, whose associated algebra is g, 
is a map J : M ^ g* such that 

di{J,X))^ixiu, (1) 

for all X G g. The previous formula reduces to Jx — ixd if there exist a 
Liouville one-form, i.e lo — dO. We also need the equivariance assumption 

{kdli,X) = {^,kdg-^X), (2) 

where Ad* denotes the coadjoint action of G in g*. 

Let us consider now an integrable system, whose configuration space is M . 
Then, there exist 2A^ functions in involution /i, I2N'- 

= 0, = 1,...,2A^, (3) 

whose differentials dij are independent at each point of an open dense subman- 
ifold K of M . The reduction procedure can be sketched as follows. Due to the 
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isomorphism between the algebra of smooth functions on M and that of the 
associated Hamiltonian vector fields, the vector fields Xi, X2N correspond- 
ing to the set of integrals in involution are in involution as well: = 0. 
Then, we can construct their symplectic action G = E?^ on the manifold. Let 
J be the associated momentum map. The theorem of Marden and Weinstein 
ensures us that, if /i is a regular value of J, then := J~^/G is a symplectic 
manifold, whose dimension is equal to dim M — N = N. 

The same construction can be used to generate new superintegrable system 
as a reduction of a given superintegrable one. In the following sections, we apply 
this approach, very general, to two physically interesting cases, like the Coulomb 
potential and the anisotropic oscillator. We point out that the integrals of the 
motion of the two reduced systems are obtained as a consequence of the reduc- 
tion scheme. The corresponding Hamiltonian vector fields are still in involution, 
since the equivariance condition ensures that the isomorphism between algebras 
described above is preserved. Also, the functional independence of the reduced 
integrals is guaranteed, if we work on a sufficiently regular manifold. In our 
specific examples, the integrals Ki, JC^v responsible for the reduction process 
are the components of the angular momentum, and the corresponding momen- 
tum mapping is J = Ki x ■ ■ ■ x Kj^. Consequently, the group action is defined 
by the cartesian product of N copies of 50(2). Finally, our reduced space is 

P^ = J-1(m)/T^, (4) 

where is the A/'-dimensional torus. As a general observation, the reduced 
Hamiltonians, as well as the integrals of motion, although keep memory of 
the original ones, are essentially different in their structure. In this sense, the 
procedure appears to be nontrivial, since it generates new dynamics from known 
ones, as will be transparent from the subsequent discussion. 



3 The Coulomb potential in even dimension and 
its reduction 

The Coulomb potential 

TV N 

^ i=i ^ i=i 

is a maximally superintegrable system, with first integrals which can be derived 
from the angular momentum and the Laplace-Runge-Lenz tensors: 

N 

Lij = XiPj - XjPi, Ai = ^pjLij - 7^, i,j = 1,...,N (6) 

i=i 

{H^,L,,}=0, {H^,A,} = (7) 

Although the whole procedure can be established in arbitrary even dimensions, 
for the sake of simplicity we will study in this work the six dimensional case, 
since the reduced system lives in a physical three dimensional space. 
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Due to the rotational symmetry of the problem, we can introduce the set of 
multipolar coordinates 



xi = xi cosaii, a;2 = a;isma;i 
2^3 = 2^2 cos 0:2 , i;4 = a:2sina;2 



2^5 = a;3 cos 0:3 , XQ — x^sinxs (8) 



and write the Hamiltonian as 



(9) 

Since the variables X4, X5 and Xe are cyclic, the corresponding momenta are 
constants and then we obtain the reduced Hamiltonian: 



ffC 1 2 7 fci ,^2 fc3 ,,^x 
3 =2^ - r + Z2+Z2+Z2 (10) 



The Hamiltonian is maximally superintegrable, that is we can find a 



This Hamiltonian lives in the reduced manifold (see Section 2), where 
the group G is the product of three copies of 50(2). 

^3^ 

set of five functionally independent first integrals (although not all of them 
in involution). Four of them can be easily obtained, being quadratic in the 
momenta [S]. However, the fifth one is not quadratic and was recently 
computed in [1^] . The method used in that work was inspired by the case when 
one of the barriers (fci) vanishes. In this situation (for instance ^3 = 0), the 
quantity (L is the usual angular momentum in the three dimensional space): 

P2Ll-piL2-2x3(-^ + ^ + |) (11) 

is an integral. The argument in [19j is that the fifth integral (which is not 
quadratic in the momenta) should reduce to this one in the limit kg, — > 0. 
This observation allows to calculate the integral we need to prove maximal 
superintegrability (see [TH] for the details) 

We will follow here a different approach, based on the theory of the reduction 
of superintegrable systems described above. The set of first integrals of the first 
Hamiltonian can be written in the new coordinates Xi. However, not all 
of them will be independent of Xi,X5,Xe, so they do not yield first integrals of 
the reduced system. Our approach consists in finding which of these quantities 
will satisfy this constraint. Equivalently, we will compute first order integrals of 
the Hamiltonian which Poisson commute with the generators of the rotation 
group L12, L34 and L^q (which are the generators of the subgroup of the rotation 
group which is used to construct the symmetry reduction). In fact, if we call 

L ~ i|rcduccdi ^ — ^|reduced (12) 

the following reduced operators are constants 

Li2 = V^i, Lm = v/2fc^, ^56 = (13) 
In fact, as we have explained in Section 2, the reduced phase space is 
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We will define a set of reduced quantities, which do not depend on X4,X5, xq 
and then, are first order integrals of the reduced system: 



E 


= H 








h 


= ^13 " 








h 


= L\,- 








h 


= -^35 " 


^ -^36 " 







(14) 
(15) 
(16) 
(17) 

{H^,h}^Q, z = 1,2,3 (18) 

We need a fifth quantity which is not quadratic. We can construct it by using 
the Laplace-Runge-Lenz tensor (remark that these quantities commute with 
ii2, ^34 and ise): 

Ti = Al+ Al (19) 
T2= Al+Al (20) 
r3= Al+Al (21) 

{i/,T,} = 0, z = l,2,3 (22) 

The seven first integrals cannot be functionally independent. In fact we only 
need a fourth order first integral. For instance, we can choose 

T = Ti + T2 + Ta (23) 

which can be written as: 



T=[-{2p'r'-ipfr) + 2r[-± + ^ + ^ + ^ +^^-^ (24) 
\Zr \ r x^ z'^ ) ) r Ar'^ 

(see [H]). We have constructed a set of five functionally independent first 
integrals (including the Hamiltonian) showing the maximal superintegrability of 
the system. The generalization of the above discussion to any (even) dimension 
is straightforward. 



4 The anisotropic harmonic oscillator in even 
dimension and its reduction 

As in the previous section, our results can be easily extended to any even di- 
mension, but we shall restrict this work to the four dimensional case and its 
reduction, the 2-dimensional case. Here we review the results contained in [T5] . 

The anisotropic oscillator in the two-dimensional case, both in classical and 
quantum mechanics, was discussed by Jauch and Hill iD], [1], [1]. The system 

N N 

1=1 1=1 

is also known to be superintegrable (in any dimension), if the ratios of the 
frequencies are rational. Let us consider a TV = 4 dimensional space and assume 

ni Til 712 "-2 
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that is, consider the Hamiltonian 



Ht = ltp' + ^Ut-'+-lt (27) 



2 2lu^ , 

i—1 i—1 i—3 

Following [IQj . we define a set of (complex) invariants, with coordinates Zi^Zi^ 
z = l,...,4: 

Zj = pj — irijUJXj, Zj = pj + irijUJXj. (28) 
It is easily checked that the expressions 

c,k = z^-zl^ (29) 

provide integrals of motion. In particular, among these integrals we have the 
angular momenta 

Lik = x^pk - XkPi (30) 

(when Ui = Uk) and the tensor 

Tik = PiPk + niUkUj'^XiXk (31) 

We will now study one of the several possible reductions of the anisotropic 
oscillator ([27]) and establish the superintegrability of the corresponding reduced 
system, computing, in an explicit way, three functionally independent constants 
of motion. 

According to the scheme adopted for the Coulomb case, let us consider the 
following change of coordinates: 

xi = xi C0SX3, £2 = xi sinxs 

X3 = a;2 cos 2:4, X4 = a;2 sina;4. (32) 
The corresponding momenta read 

P1=-P3^^+P1C0SX3, P2=P3^^^+PlSmX3, 
P3^-P4'-^+P2COSX4, P4^P4^^ +P2SmX4. (33) 

and the system in the new coordinates is: 



Ht -UpI+pI) + ^ + A + ^^'{nixi + nixi) (34) 



Setting 

P3 = V^) P4 = \/fc2 (35) 

we get the reduced system: 

H2 = 2(^1 +^2) + ^ + ^ + [n^x^+n^x^ j (36) 

Notice the appearance of Rosochatius-type terms [TB] in the reduced system. 
As in he Coulomb case, this Hamiltonian is defined in a reduced manifold 
and the group used in this reduction is 50(2) x 5*0(2). 
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The integrals of the system (P7|. given in will be used to construct a 
set of integrals for the reduced system. Those that will survive the reduction are 
the ones that are left invariant by the 50(2) x 5*0(2) rotations. The arguments 
are the same we used in the Coulomb potential, in fact, they must Poisson 
commute with 

■^12= 2;rnj(-^i^2 - 22^i) = iiP2 - ^2^1, 

■^34= 2K^(^3Z4- Z4Z3) = X3P4- XiPs. (37) 

The computation of these commuting quantities is not as straightforward as 
in the Coulomb potential, where we found the integrals by simple inspection. 
We will use z coordinates to solve our problem. The Poisson bracket can be 
written in terms of the Zi variables as 



N 2k 



{/(z„z-,),.9(z.,z-.)} = -2ia;5: ^ (38) 

k=l j=2fe-l \ 3 3 J 3 



hence, functions of Zk, Zk Poisson commuting with L12 and L34 must satisfy 

Z2dzJ - Zldz^f + 22921 / ~ Zldg^J = 0, 

Z4dz3f - zsdzj + Z4dzj - zsdg^f = 0. (39) 
A basis for the corresponding 50(2) x 50(2) invariants is given by 

^l=Zl+z|, ^i=Zi+z|, TJl = ZlZl + Z2Z2, 

^3^zl+zl, (s^zl + zj, 7]2 ^ Z3Z3 + Z4Z4, (40) 

These rotationally invariants will be constants of motion of our system if they 
satisfy: 

{Ht f i^i, Cum, ^3, m)} = 0. (41) 

The following expressions are solutions of equation (|¥T|) 

E^ = ^i\z^\' + \z2n E2 = \i\z,\^ + \Z4\'), 

Ql = (zf + zlY- [zl + Z|)"i , Oi = (f 2 + (^2 ^ ^2)„, ^ (42) 

h = [zl + zl){zl + 4), I2 = {zl + zl){zl + zl). 

however, only five of these integrals are functionally independent. 

We now perform the reduction for the integrals of motion using the change 
of variables ([32]) . (l33|) . The integrals ((ST]) reduce to constants L12 — -y/fci, -^^34 — 
\fk2. The first four integrals in (|42p reduce to nontrivial integrals for the Hamil- 
tonian p6p . namely 



Ex 


1 2 


^ 2x1 


1 1 2 2 2 




E2 


1 2 

= 2^2- 


^—2 
2x2 


1 1 2 2 2 
H — n^LO Xn 

2 




Qi 




ki 
'xl 


n\uj^x\ — 2iniijjpiXi)'^^ {p2 H 


k2 

X2 


Qi 




ki 
' xl 


n\uP'x\ + 2iniijjpiXi)"'^ {p2 H 


_k2 

x\ 



nlu'^xl + 2in2wp2a;2PS,) 
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The remaining two integrals in (j42p give nothing new 

h - ^{El - kmluj^), h = 4(£;| - k2nlLu^). (44) 

Three functionally independent real integrals of motion of the system with 
Hamiltonian (|27|1 can be chosen to be 

{E^,E2M^\{Qi+Qi)}- (45) 

They are of order 2, 2 and 2(ni + n2) in the momenta, respectively. Their 
existence is the proof of the maximal superintegrability of the considered system. 

The integral of motion Q simplifies to give a second order one in two cases 
(that were known previously [7], [10]). They are 



I) 111 — n2 — I 



4:EiE2 - Q , ,2 , kixl k2xl 

^{piX2-p2Xl) +-^ + ^- (46) 



II) ni = 1, 712 = 2, fca = 
(%ElE2-Q 



kiE, 



1/2 



Pl{x2Pl - X1P2) - Uj'^x\x2 + k 



X2 



(47) 



The integrals (|46|) and (|47p are responsible for the separation of variables in polar 
and parabolic coordinates, respectively. The integrals {£'i,i?2} are responsible 
for the separation in cartesian coordinates. 



5 Conclusions 

From the previous examples it is clear that it would be very interesting to con- 
struct, in a systematic way, transformations mapping a superintegrable system 
into another system, that is also superintegrable, and defined in a reduced phase 
space. The role of higher order groups of transformations generated by the flow 
associated to integrals that are polynomials in the momenta remains to be fully 
investigated. 

We also mention that an algebraic approach to superintegrable systems, in- 
dependent of the geometry of the phase space, and based on a suitable nilpotency 
condition, is presently under investigation fS]. 

A quantum mechanical version of this reduction procedure is also to be 
understood. 
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